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The  Nonlinear  Dynamic  Response  of  an  Elastic-Plastic  Thin  Plate 
Under  Impulsive  Loading 


by  Uang  Xintian,  Hong  Shantao,  and  Weng  Zhiyuan 
Tongji  University 


Paper  received  on  August  25,  1983. 


n  this  paper  the  effects  of  the  physical 
and  geometrical  nonlinearities  in  a  thin  plate 
are  treated  as  equivalent  body  forces  and 
equivalent  loads.  Using  the  concept  of 
influence  functions,  an  analytical  method  for 
the  thin  plate  problem  with  both  kinds  of 
nonlinear  effects  is  presented.  In  theoretical 
analysis,  the  influence  of  plastic  region  which 
extends  in  depth  is  considered  and  an 
incremental  formula  of  plastic  strain  is  derived 
by  applying  the  incremental  plastic  theory.  In 
the  calculation  of  practical  examples,  the 
numerical  solutions  for  nonlinear  dynamic 
responses  of  an  elastic-plastic  thin  plate  are 
obtained  for  various  hardening  coefficients  and 
different  impulsive  loads,  all  of  the  results 
are  quite  regular.  (C Tol o 
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I .  Preface 


Due  to  the  requirement  of  the  antiknock  and  anti-seismic 
design  of  a  structure,  the  nonlinear  dynamic  response  analysis  of 
a  structure  has  been  widely  recognized. 

Two  different  kinds  of  nonlinear  problems  will  be  encountered 
in  the  structure  dynamic  analysis  if  the  structure  undertakes 
short-duration  and  large  impulsive  loads.  One  is  the  physical 
nonlinearity  which  is  the  nonlinearity  of  the  relationship  between 
material  stress  and  strain.  The  other  one  is  geometrical 
nonlinearity  which  is  the  nonlinearity  of  the  relationship  between 
strain  and  displacement  which  occurs  while  the  structure 
deflection  exceeds  the  small  deflection  range.  In  the  engineering 
application,  more  attention  is  paid  to  the  dynamic  response 
analyses  which  cover  both  kinds  of  nonlinearity.  It  is  a  rather 
complicated  problem.  Generally,  using  numerical  computation 
methods  is  the  only  way  to  obtain  the  approximate  solution.  Among 
those  numerical  methods,  the  finite  element  method  is  the  most 
effective  one.  The  method  presented  in  this  paper  is  an 
"influence  function"  numerical  method.  It  is  conceived  from  A.  A. 

JC7  *yushchin ’s  elastic  method  for  a  plastic  problem,  A.  5. 

M»/*mir’s  method  for  a  large  deflection  problem,  and  T.  HI. 

Lin's  method  for  a  dynamic  problem.  Furthermore,  the  concept  of 
influence  functions  is  applied  to  solve  the  dynamic  response  of  a 
thin  wall  problem  with  both  kinds  of  nonlinear  effects.  Compared 
with  the  finite  element  method,  the  method  presented  in  this  paper 
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is  simpler.  Its  computational  time  is  much  less  than  the  finite 
element  method.  Its  precision  is  higher  also. 

II.  Basic  Assumption  and  the  Influence  Function  of  a  Simple 
Suppo rted  Plate. 

1 .  Basic  Assumption 


1)  The  loading  is  a  suddenly  increased  uniform  load 
which  varies  with  time. 

2)  The  structure  is  an  elastic -plastic  thin  plate.  Both 
elastic  and  plastic  regions  follow  the  Kirchhoff 
assumption . 

3)  The  material  of  the  structure  contains  the  dual-linear 
hardening  characteristic  as  shown  in  Fig.  1. 

4)  The  plastic  theoretical  analysis  adopts  the  isotropic 
hardening  model.  The  material  follows  the  von  Mises 
yielding  function  and  the  associated  Prandtl-Reuss 
flow  rule. 

5)  The  inertia  force  effects  at  both  x  and  y  directions 
are  neglected  in  the  dynamic  response  calculation. 
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2.  The  Influence  function  of  a  Simple  Supported  Plate 


The  motion  differential  equation  of  a  thin  plate  can  be 
written  as 

+  _£*_*«,(*)  <1 

9 

Its  resonance  frequency  is 


Assuming  its  particular  solution  is 

■••EE  Xm(x )  Y.(y)  r.(f)  ( 3 } 

■•I  •-! 

After  normalization,  choose  Xm(x)  - *?*-,  Y.  (y)  lin  — 

Applying  the  orthogonality  of  the  trigonometric  function,  and 
choosing  the  initial  condition  as  rnmum0  when  t=0,  we  have 


r-(0 >-  <4> 


The  dynamic  response  of  a  simple  supported  rectangular  plate 
due  to  the  point  load  of  an  unit  strength  at  (£,  9)  is 


mix^y 


•"■-a-SS - 1 - te— * -  (s) 


X  <1 


This  is  the  Green  function  of  the  dynamic  response  of  a  simple 
supported  thin  plate.  After  the  time  and  special  variables  in  the 
Equation  5  are  separated,  the  dynamic  influence  coefficient  of  a 
thin  plate  deflection  is 
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x  Cl  -  cos <*>«.(/»•  AO D 


(6) 


Equation  6  represents  the  deflection  response  at  grid  (i,  j)  when 
time  step  is  equal  to  p  due  to  the  step  load  of  a  unit  strength 
act  on  grid  (k,  1)  when  t=0. 


Similarly,  the  static  influence  coefficient  of  a  thin  plate 
deflection  can  be  obtained  from  the  thin  plate  static  differential 
equation . 

S(i,  £ 

sin  -E-«Ay> 

Both  dynamic  influence  coefficient  6  and  static  influence 
coefficient  S  of  a  thin  plate  are  called  influence  function  for 
short  in  this  paper. 


sin-52-  (i  Ax)  sin-?r-(/Ay)  sin-55,  (k&x) 

0  U  9 

(”r 


III.  The  Dynamic  Responses  of  an  Elastic-Plastic  Thin  Plate  Small 
Deflection  When  Hardening 


The  non-hardening  case  of  the  dynamic  responses  of  an  ideal 
elastic-plastic  thin  plate  small  deflection  has  been  discussed  in 
Ref  (3).  This  section  will  depict  the  dynamic  response  of  an 
elastic-plastic  thin  plate  small  deflection  while  considering  the 
hardening. 


S 


The  motion  equation  of  an  elastic-plastic  thin  plate  (small 


deflection)  is  written  as 


OV4m>  + 

0 


(8) 


where  q(x,y,t)  is  the  lateral  equivalent  loading  caused  by  the 
plastic  strain  inside  the  plate,  and  is  written  as 


<«i+  + vt[)idz 

Oxdy  \e’’*dz\ 


+  2(1  -*)- 


(9) 


where  el,  ej,  ei  are  time  varied  plastic  strain  components. 


The  thin  plate  is  divided  into  uniform  grids  while  the  time 
history  is  divided  into  equal  intervals.  Considering  the  case  of 
simple  supports  of  four  sides  of  the  plate,  substitute  the  dynamic 
influence  coefficient  G  shown  in  Equation  6  into  it,  then  the 
dynamic  response  successive  substitutional  equation  of  a  thin 
plate  at  time  pAt  is 

«  H 

2  £*»(••/**♦*.£>  CAq(k,l,o)  +  Aq(,k,  /,o)3 

i.i  i*i 

♦  I]  H  L.S  ti,i,k,l,p-r)CAq(k,lfr)  +A9(*,/,r>D 

f-i  t.i  i.i  (10 

+  £  'jrt—G(i,i,k »M)  CAq<** I »P>  +  Ay(M./>>3 

Where  Aq(kfl,r)  is  the  lateral  loading  increment  act  on  grid  (k, 
1)  at  ,  and  Aq{k,l,r)  is  the  plastic  equivalent  loading 

increment  act  on  grid  (k,  1)  at  while  hardening  is 

considered.  G  (i,jjkj,p-r)  --jC(?(i,  j,k,l,p-r)  +  GU,j,h, /,*-r+  1)3. 
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Considering  the  dual -linear  situation  of  the  material 
stress-strain  relations  (see  Fig.  1),  a  is  a  hardening 
coefficient,  a,  is  a  plastic  yielding  limit.  Using  the 
Prandtl-Reuss  flow  rule  and  the  broad  sense  of  the  Hooke  law,  the 
relationship  between  pre-plastic  ized  strain  increments  Ael,,  Ae'„.  Ae'„ 
and  full  strain  increments  At.,,  Ae„,  Atm  at  time  pAt  is 


Ae!,*-y-C(^  +  rB>Ae.,+  (B  +  rA)Ae„  +  (1  -  *)CAe,„] 
A  *;,  «  -y-c  (A  +  rB)  At.,  +  (B  +  pA)  Ae„+ (1-p)  CAt.^J 
A*'.„  -  -§rC  (A  +  vB)  At.,  +  <fl  +  vA)  At,,  +  (1  -  *)  CAe„P 


(11) 


where 

C*3tw-,  (12) 

r  -  (1  -  +  A(A  +  „B)  +  B(B  +  vA)  +  -irz!  c» 

2 


the  stress  strength  o,p_l  is 


o„-x  -  <<*!,-■  +  a!,.,  +  c„.,awM  +  3rJw-,)  *  (13) 

The  increment  form  of  the  plastic  strain  equivalent  loading  q  is 


Aq  ■  -■y-~“r£*g“rJ<A«I  +  *A«J)«fe +-j£jr^{A*\  +  vAt'.)xdg 

m  t  *  (14) 

Notice  that  the  plastic  strain  increment  formula,  Equation  11, 
derived  in  this  paper  which  considers  the  hardening  situation,  is 
consistent  with  the  matrix  forms  used  in  the  elastic-plastic 
finite  element  methods  which  are  presented  in  Ref  (161  ,  (181  . 
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The  form  of  the  dividing  point  £  ( x,y,t )  between  elastic  and 


plastic  along  the  depth  in  the  plate  can  be  expressed  as 


€  (*,y,  t) 


.  (  E  r,Vw  Vw,*  ,d*u»  .  &W,1  ,&w_  . 

+  +  ‘r*r)  (e^ 

/  ym  .  &w  \  .  .  ,  0*111  v  * lil  -* 

^+vair)  +  3(1-*')  (-557)  JT} 


0*u>  \ 

’ra7) 


(  15) 


During  the  successive  substitution  process,  based  on  the 
Kirchhoff  hypothesis,  the  full  strain  increment  can  be  obtained  by 
using  the  formulas  below 


A  if,,*  -x- 
{  A«r„  *  -  * 

A  if,,,  ~ 


0*(Au>,) 
0x* 

_  9*  (Aw,) 
9t/1 

_9*(Auy) 

dxdy 


(.  16) 


where  Aa,  =  a,-»,^.  The  total  stress  at  time  p-At  is 


<7,,  *  Q af- 1  ^  A<7„ 

1 +  Act„ 

Ar„, 


(17) 


For  the  convenience  of  comparison,  the  same  rectangular  plate 
used  in  Ref  131  is  selected  to  calculate  the  dynamic  response  of 
an  elastic-plastic  thin  plate  small  deflection  when  the  hardening 
occurs . 
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The  thin  plate  is  shown  in  Fig.  2,  where  a=25.4cm,  h=1.27cm, 

5  2  3  2 

E=7. 03x10  kg/cm  ,  *>=0.3,  r  =  2.77  g/cm  ,  =2109kg/cm  , 

2 

and  its  four  sides  are  simple  supported.  q=21kg/crn  ,  the  time 

interval  is  chosen  as  1/48  of  the  basic  period  T  , 

-4  1 

At=0. 223x10  8.  Take  the  advantage  of  symmetry,  divide  the  one 

quarter  of  the  plate  into  4x4  grids,  and  divide  five  intervals  in 

the  depth  of  the  plate.  Based  on  four  different  conditions  <*=  o, 

0.2,  0.5,  and  1.0,  the  results  of  calculation  are  plotted  in  Fig. 

3.  From  Fig.  3,  it  shows  that  the  maximum  deflection  decreases 

while  a  increases,  and  the  thin  plate  appeared  to  have  a  hardening 

tendency.  The  curves  are  quite  regular. 


<■  r 

•i 

r 

l 

1 

«# 

_ 1 _ 

« 

» 

Fig.  2 


Fig.  3 


When  the  load  is 
the  small  deflection 
nonlinearity  effects 
will  be  discussed  in 


further  increased,  the  deflection  will  exceed 
range.  In  this  case,  the  geometrical 
must  be  considered  at  the  same  time.  This 
the  next  section. 
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IV.  Analysis  of  a  Problem  which  Contains  both  kinds  of 

Non  1  inear i ties 

In  the  anti-knock  and  anti-seismic  design  of  a  structure,  the 
impulsive  loads  are  usually  violent.  This  often  introduces  both 
physical  and  geometrical  nonlinear  effects  in  the  structure. 
Therefore,  this  sort  of  problem  has  obvious  engineering 
significance. 

In  order  to  systematically  analyse  the  problem,  we  shall 
attack  the  topic  start  from  the  static  part  then  extend  to  the 
dynamic  part. 


1.  The  static  analysis  of  an  elastic-plastic  thin  plate  large 
def lection 


The  coordinate  system  is  shown  in  Fig.  4.  Based  on  the 
Kirchhoff  hypothesis,  we  have 


du  1  /  9u> 

dx  2  {  dx  ' 


__  9*tt> 

e*‘  =  e"  ~  '  dx 1 


dv  1  /  dw  t Vw  d*w 

*'*  dy  +  2  '  dy  >  ~  *W  dyz 


1  ,  du  dv  dvu  dw  v  _  ym 

*•»*  2  ^  dy  *  dx  dy  dx  ’  dxdy 


€‘*~  *  dxdy 


(  18) 


where  e,t%  are  plane  strains.  The  total  strain  is  composed  of 
the  elastic  strain  and  the  plastic  strain.  For  the  plane  stress 
problem,  the  stress-strain  relation  is 
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Fig.  4 


Substitute  Equation  18  into  Equation  19,  and  then  integrate  it 
along  the  entire  plate  depth,  and  we  have 


1u  .  dv  J/  9u’\*  .  v  /  dtv  v'1  E 

2'-dr>  +T(v)  3  1-*-* 

J)  dt  =  N.,-N, 

)ti  An  1  Aim  v  *  it  ,  Atm  ,  h  /** 


where  JVH,  AT,,*  are  the  film  forces  of  the  corresponding  middle 
plane  strains,  while  77,,  27„are  the  film  forces  of  the 
corresponding  plastic  strain  in  the  large  deflection  condition. 
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The  static  equation  of  an  elastic-plastic  large  deflection  is 


r0V4w  *  9  +  $  +  ff1 

m  ~fi  +  f1  *■  22 ) 

where  q  is  the  lateral  load,  q  is  the  plastic  strain  equivalent 
1 

load,  and  q  is  the  large  deflection  lateral  equivalent  load. 


_i_  d1^  6*u>  .  0*4  0*w  d*m 

’  6yl  Ox *  Ox*  0yl  bxOy  bxOy 


<23) 


F  is  the  plastic  strain  equivalent  body  force  of  the  stress 
function,  written  as 


T  ■  -  ♦  -&1  '<*  - 1  ]  <“> 

1 

F  is  the  large  deflection  equivalent  body  force  of  the  stress 
function,  written  as 


A 

The  relationship  between  the  stress  function  ^  ,  the  middle  plane 
film  force,  and  the  middle  plane  strain  is 


N* 


*.r 


dV_ 

bxdy 


(26) 


.  _  a** 
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The  successive  substitution  formula  can  be  written  as 

wAi,  />«  E  £ S /,*,/)  C Aq,  (k,l)+  Aq! (*./> +Aq,<*./> 3  <28) 

»-i  i-i 

Solutions  can  be  obtained  by  solving  the  Equation  (28)  and  the 
second  equation  of  Equation  (22)  simultaneously.  The  entire 
process  can  be  stated  briefly  as  follows:  The  load  q  is  increased 
by  the  same  increment  Aq,  the  rth  load  increment  is  called  the 
rth  load.  Let  Aq!  =  A?,  =  AF,  =  &F1,  =  0  ,  substitute  it  into  Equation  (28) 
and  obtain  w?  .  This  u>?  can  be  used  as  the  initial  value  of 
successive  substitution.  The  first  approximate  value  of  the 
plastic  strain  increment  of  this  stage  can  then  be  calculated  by 
considering  the  total  stress  created  under  the  final  load  of 
previous  stage  load,  and  the  load  increment  Aq  of  this  stage. 
FindAF!-  AFt'and  substitute  them  into  the  second  equation  of 
Equation  (22),  then  A$.  can  be  obtained  by  applying  the  finite 
difference  method,  thus  Aq!'  and  Aq?  can  be  calculated.  Again 
substitute  both  values  into  Equation  (28)  to  obtain  the 
approximate  value  after  the  first  successive  substitution  *»?.... 
Repeat  this  process  until  the  converge  cr i ter iaAwJ-AwV1! <eis  met, 
which  means  the  difference  is  negligible.  Then  add  one  more  load 
increment  Aq.  Repeat  the  entire  above  processes  until  the  final 
load  reaches  q.  The  final  solution  w( i ,  j)  is  the  static  solution 
of  the  elastic-plastic  thin  wall  large  deflection  under  the  load 

q* 
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This  paper  takes  the  square  plate  shown  in  Fig.  2  as  a 
practical  example  and  conducts  the  calculation  by  treating  it  as  a 
static  problem.  The  boundary  condition  is  that  the  four  sides  of 
the  plate  are  hinge-supported.  This  is  the  simpler  case.  Four 
difference  conditions  are  considered  in  the  calculation:  Elastic 
small  deflection  (E.S.S),  Elastic  large  deflection  (E.L.S), 
Elastic-Plastic  small  deflection  (E.P.S.S),  and  Elastic-Plastic 
large  deflection  (E.P.L.S).  Four  resulting  curves  are  plotted  in 
Fig.  5. 


For  the  case  of  elastic  large  deflection  problem  of  a  thin 
plate,  A.  C.  Vol"*ie  gave  the  first  stage  approximate 
calculation  formula  of  a  square  plate, 

*•♦*-€* 

Under  the  four-side  hinge -supported  condition,  A=7 .5 ,  B=22.  The 
nondimensional  parameter  ***!?.? ***  *  where  f  is  the 

maximum  deflection.  The  more  precise  solution  should  be  located 
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at  the  left  side  of  this  approximate  solution,  Ftef  (41  . 
Meanwhile,  based  on  Ref  1101,  this  paper  obtains  one  curve  by 
using  interpolation  method.  From  Fig.  6,  it  shows  the  method 
presented  in  this  paper  has  higher  precision. 


Key:  (1)  Method  in  this  paper;  (2)  Ref  [41; 

<3)  More  precise  solution;  (4)  Ref  [101 

Fig.  7  and  Fig.  8  show  the  plastic  region  distribution  on  the 
upper,  lower  surface  of  the  plate  and  along  the  depth  of  the  plate 
of  the  case  which  considers  large  deflection  effects.  Fig.  9 
shows  the  deflection  curves  at  the  different  cross  sections  of  the 
plate  which  again  are  the  results  of  the  four  different  cases. 
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2.  The  dynamic  response  analysis  of  the  elastic-plastic  large 
deflection  of  a  thin  plate 


The  motion  equation  of  the  elastic-plastic  large  deflection  of 
a  thin  plate  can  be  written  as 


9  (31  > 

v>  -  F  ♦  F 


Compared  with  the  static  equation,  the  above  equation  contains 

only  one  term,  an  inertia  term,  more  than  Equation  (22).  Other 

than  the  time  factor  need  to  be  added,  the  expression  forms  of  the 

stress,  strain,  internal  force,  internal  torque  as  well  as  q, 

1  —  1 

q  ,  F,  F  ,  +  and  etc.  for  the  thin  plate  in  the  dynamic  case 
are  similar  to  those  in  the  static  case.  In  this  paper,  the 
inertia  term  is  handled  by  applying  T.  H.  Lin’s  dynamic  problem 
solving  method  and  introducing  the  dynamic  influence  coefficient 
G(  i ,  j,  k ,  1,  p)  of  a  thin  plate  deflection.  The  successive 
substitution  formula  is 


•  2D  £  C(i,  j,  ft,  i»#)CA«  (*,(,•>  <*’!.»)  +  6c'(i,1.0)) 

••ft  l-l 

f  -  •  «  ft 

♦  E  E  E ^ /* KUp-r) CAf (ft, /. r)  ♦  64 ♦  Ag'(A,(,r>] 

*•1  1*1  I •  I 

tt  ft 

♦  E  E  j-C(i,;.ft.M>CAf(ft,I,^>  ♦*<*,/.*>  ♦  ^<*.1.^)1 

»•!  1*1  * 


(32) 


where  gf'UM.r)  is  the  large  deflection  equivalent  load  increment 
which  acts  on  grid  (k,  J)  at  time  #»#,  .  The  finite  difference 
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incremental  form  of  the  second  equation  of  Equation  (31)  on  the 
grid  <i,  j)  at  time  pAt  is 

*  +  ^  +  4^1-11/  (33) 

+  -  »)  *  if  l«tl  *  if  fill 

The  successive  substitution  process  is  similar  to  the  previous 
case.  For  the  convenience  of  comparison,  the  square  plate  shown 
in  Fig.  2  is  still  used  as  a  practical  example  to  calculate  its 
dynamic  response.  The  boundary  is  still  the  hinge -supported 
points . 

The  impulsive  load  is  a  suddenly  increased  platform  load. 

2 

Three  difference  loads  which  are  q=28,  35,  56  kg/cm  , 

respectively,  are  used  in  the  calculations.  Four  different  cases 

are  analyzed  in  this  paper;  they  are:  elastic  small  deflection 

(E.S.D),  elastic  large  deflection  (E.L.D),  elastic-plastic  small 

deflection  (E.P.S.D)  and  elasti c -plasti c  large  deflection 

(E.P.L.D).  The  dynamic  response  curves  are  then  obtained  as  shown 

in  Fig.  10,  11,  and  12.  In  order  to  discuss  the  geometrical 

nonlinear  effect,  the  impulsive  load  was  increased  to  q=70 
3 

kg/cm  in  this  paper.  The  cases  of  E.S.D,  E.L.D,  and  E.P.L.D 
were  calculated.  The  dynamic  response  curves  for  these  three 
cases  are  shown  in  Fig.  13. 

This  paper  also  calculates  the  cases  when  the  impulsive  load 
is  a  triangular  and  a  square  function.  The  associated  dynamic 
response  curves  are  shown  in  Fig.  14  and  15. 
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From  these  curves,  we  notice  that  the  thin  plate  is  "softened" 
when  the  plastic  effect  is  considered,  while  the  thin  plate  is 
"hardened"  when  the  geometrical  nonlinear  effect  is  considered. 
Furthermore,  the  geometrical  effect  abruptly  increases  along  with 
the  increasing  of  the  impulsive  load.  The  curves  show  the 
particular  regularity. 

Due  to  there  being  no  literature  which  discusses  the  influence 
of  plastic  region  which  extends  in  depth  and  also  considers  the 
thin  plate  dynamic  response  with  both  kinds  of  nonlinearity 
effects,  a  direct  comparison  and  discussion  can  not  be  made. 
However,  a  "degraded"  situation  is  made  in  this  paper  to  conduct 
an  indirect  comparison  and  discussion. 

Ref  1111  applies  the  finite  difference  method  to  the  studying 
of  the  dynamic  response  of  a  thin  plate  elastic-plastic  small 
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T 


deflection  problem.  That  paper  uses  nondimension  parameters.  For 
the  case  of  the  four-side  simple  supported  square  plate,  we  have 


P-'T. TOST  T'  1,73 

?  r _ 6  „ _ t 


PS 


2  nal 


.  m 

nr 


<341 


’  p,  °  qi,  '  r, 

where  M ,  -  is  the  thin  plate  yielding  moment.  Notice  that  Ref 

[111  uses  the  finite  element  method  to  obtain  T  which  is 

1 


inconsistent  with  the  result  obtained  by  the  analytical  method 
presented  in  this  paper.  The  period  used  in  Ref  [111  was 
converted  into  the  analytical  solution  value  when  plotted  in  Fig. 

t 

16.  For  the  square  plate  shown  in  Fig.  2,  q  =2  is  equivalent  to 

2 

q=54.8kg/cm  .  From  Fig.  16,  the  curve  obtained  by  the  method  in 
this  paper  in  terms  of  calculating  the  case  which  degraded  to  the 
elastic-plastic  small  deflection,  matches  quite  well  with  that, 
obtained  by  using  the  finite  element  method  presented  in  Ref 
[111 .  Both  curves  appear  even  closer  for  the  maximum  deflection. 


i 


Fig.  16 
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If  the  plastic  effect  is  neglected  in  the  case  of  two 
nonlinearities,  the  problem  will  "degrade"  to  an  elastic  large 
deflection  condition.  The  results  of  such  a  "degraded"  case  by 
using  the  method  presented  in  this  paper,  therefore,  cart  compare 
with  those  presented  in  Ref  [61  . 


Let's  introduce  the  related  parameters  in  Ref  £61 , 


(35) 


Ref  £61  points  out  that  the  maximum  amplitude  A  of  the 

max 

dynamic  response  of  an  elastic  thin  plate  large  deflection  in  fact 
is  one  of  the  real  roots  of  Equation  E(A)  *  A*  + 

Obviously,  the  solution  based  on  the  elastic  small  deflection 

3 

theory  can  be  obtained  if  A  =0.  This  paper  uses  to 

represent  the  maximum  deflection  of  an  thin  plate  elastic  dynamic 
response  which  is  obtained  according  to  the  large  deflection 
theory  and  small  deflection  theory,  respectively.  It  must  be 
pointed  out  here  that  only  the  first  term  (N=l)  of  the 
trigonometric  series  is  taken  in  Ref  £61  which  implies  that  the 
thin  plate  is  analyzed  by  treating  it  essentially  as  a  single 
degree  of  freedom  system.  Both  u>su  are  the  first  stage 

approximation  only.  This  paper  applies  the  method  of  influence 
function  where  N=7  is  taken  in  the  trigonometric  series.  The 


results  are  listed  in  Table  1. 
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Table  1 


1>«  * 

(  k*  ) 

'  cm2  f 

w$  n  (cm) 

w  Lit  (cm) 

ntiin 

■ 

yflfl 

(4) 

35  ! 

i  0.923 

; 

0.881 

4.85 

|  0.845 

0.813 

i 

5.06 

42 

1.108 

|  1.056 

4.88 

1.001 

0.953 

5.09 

56 

1.477 

1.408 

4.89 

1.264 

1.205 

4.87 

70 

1.846 

1.761 

4.88 

1.494 

■a 

4.68 

Key:  (1)  Load;  (2)  Ref  161 ;  (3)  Method  in  this  paper; 

(4)  Error  . 


In  order  to  further  eliminate  the  effect  of  the  grid  division 
in  the  finite  difference  form  and  the  number  of  the  term  in  the 
trigonometric  series,  this  paper  introduces  a  large  deflection 
influence  coefficient  K.  : 


H.*U>Lii/WsM  (38) 

This  parameter  reflects  the  reducing  rate  of  the  maximum 
deflection  of  a  thin  plate  dynamic  response  when  a  large 
deflection  effect  is  considered.  Thus,  Table  2  is  obtained. 


Table  2. 


<1)«  « 

(J^r) 

'em* ' 

(2) 

I 

(3) 

(4) 

ffi  fi  ("o' 

35 

0.9252 

0.9233 

0.21 

42 

0.9037 

0.9019 

0.20 

S« 

0.8555 

0.1557 

0.02 

70 

0.0091 

0.1106 

0.18 

Key:  <1)  Load;  (2)  k-  obtained  by  Ref  161  ; 
(3):  obtained  by  this  paper;  (4)  Error. 
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From  both  Table  1  and  2,  the  dynamic  responses  of  an  elastic 
large  deflection  obtained  by  using  the  influence  method  presented 
in  this  paper  have  higher  precision. 

This  paper  also  calculates  the  case  of  choosing  N=1  in  the 
trigonometric  series,  and  compares  the  result  with  Ref  [61 .  It  is 
found  that  the  error  can  be  further  decreased  to  2-3%. 

In  the  practical  example  calculation  of  finding  the  dynamic 
responses  of  an  elastic-plastic  large  deflection,  when  even  the 

time  interval  At  is  doubled,  the  difference  of  the  maximum 

2 

deflection  will  not  exceed  10%  in  the  q=56kg/cm  case.  This 
states  that  the  convergence  of  this  paper's  method  is  better  also. 

In  a  word,  the  thin  plate  dynamic  response  curves  of  different 
conditions  obtained  by  using  the  method  presented  in  this  paper, 
contain  the  particular  regularity.  Because  no  compatible 
literature  was  found,  this  paper  compared  the  results  of  the 
"degraded"  condition  with  the  corresponding  literature.  The 
results  show  remarkably  well.  This  explains  indirectly  that  it  is 
feasible  to  analyze  the  problem  with  both  kinds  of  nonlinear 
effects  by  using  the  influence  function  concept. 
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Discussion 


V  . 


1)  This  paper  applies  the  concept  of  influence  function, 
considers  the  influence  of  the  plastic  region  which  extends  in 
depth  and  in  the  x,  y  direction  of  the  plate,  and  analyzes  the 
thin  plate  dynamic  response  problem  with  both  kinds  of  nonlinear 
effects.  Through  calculation  and  comparison  of  the  results  of  a 
degraded  case  with  the  related  literature  practical  examples,  it 
indirectly  proves  that  the  method  presented  in  this  paper  is 
feasible.  From  the  given  elastic  solution,  this  paper  solves  the 
influence  function.  For  those  more  complicated  problems,  if  their 
elastic  solutions  are  not  given,  they  can  be  found  through  finite 
element  methods  as  pointed  out  in  the  Ref  131  abstract. 

Basically,  the  method  presented  in  this  paper,  therefore,  can  be 
applied  to  the  nonlinear  dynamic  response  analysis  of  complicated 
structures  whose  elastic  solutions  are  given.  The  calculation 
time  takes  around  35  minutes  for  each  case  in  a  719  machine. 

2)  It  is  noticed  through  calculation  in  examples  that  the 
physical  nolinearity  makes  a  thin  plate  "soften"  while  the 
geometrical  nolinearity  makes  it  "harden".  Although  they  both  are 
coupling  each  other,  their  effects  on  the  structure  are  contrary. 
Along  with  increases  of  the  impulsive  loads,  both  nonlinear 
effects  build  up;  however,  the  effect  of  geometrical  nonlinearity 
builds  up  abruptly.  From  the  response  curves,  we  noticed  that  the 
majority  of  a  structure  E.P.L.D  curve  fell  inside  of  the  E.S.D 
curve  when  the  impulsive  load  reached  a  certain  level.  N.  Jones 
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pointed  out  that  the  plate  strain  rate  is  the  primary  factor  when 
the  impulsive  loads  are  small,  while  the  film  force  of  the  plate 
middle  plane  is  the  primary  one  when  the  impulsive  loads  are 
large.  However,  it  lacked  precise  numerical  results.  This  paper 
has  now  given  the  quantitative  results  and  the  associated 
conclusions.  Moreover,  they  are  consistent  with  the  viewpoint  of 
N.  Jones. 

3)  After  considering  the  geometrical  nonlinear  effect,  the 
neutral  plane  of  a  thin  plate  deviates  from  the  plate  middle 
plane.  The  plastic  strain  of  the  plate  extends  in  depth,  its 
distribution  is  not  symmetrical  with  respect  to  the  plate  middle 
plane.  The  examples  show  that  such  a  situation  can  cause  the 
reduction  of  the  elastic  core  in  the  plate  layer,  thus  forming  a 
plastic  yielding  region.  This  phenomenon  is  very  significant- 
while  studying  the  development  regularity  of  a  thin  plate  plastic 
region  or  a  plastic  hinge  line  along  with  time. 

A)  When  the  impulsive  loads  are  large,  under  a  finite 
deflection  situation,  we  notice,  from  the  calculated  curves,  the 
ratio  between  the  difference  of  two  maximum  deflections  obtained 
from  a  thin  plate  E.P.S.D  curve  and  E.S.D  curve,  and  the  maximum 
deflection  of  E.S.D  curve,  far  exceeds  the  value  under  the  30% 
conclusion  drawn  by  some  literature  which  is  based  on  the  small 
deflection  theory.  When  the  impulsive  load  is  large,  the  thin 
plate  dynamic  response  has  exceeded  the  small  deflection  region. 
Therefore,  the  geometrical  nonlinear  effect  has  to  be  taken  into 
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account;  otherwise,  the  error  will  abruptly  enlarge  along  with  the 
increase  of  load.  This  is  also  where  the  significance  of  this 
paper  lies. 

5)  The  nonlinear  dynamic  response  of  an  elastic-plastic  thin 
plate  studied  in  this  paper  belongs  to  the  cases  where  the  ratio 
of  the  maximum  deflection  and  the  thickness  of  plate  is  not  too 
large.  Because  when  the  deformation  of  an  elastic-plastic  plate 
is  large,  the  precise  analysis  becomes  complicated,  and  it  is 
beyond  the  scope  of  this  paper. 

6)  The  method  presented  in  this  paper  can  only  obtain  the 
approximate  solutions,  because  this  method  adopts  the  following 
approximations  to  find  out  the  numerical  solutions.  They  are  a 
successive  substitution  method  to  approaching  the  solution  step  by 
step,  a  time  and  special  variables  separation  method,  finite 
difference  forms  as  well  as  the  method  of  finding  a  plastic  strain 
increment  of  the  current  step  by  using  the  final  full  stress  of 
the  previous  step,  and  so  forth. 
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